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ABSTRACT 

A novel method for generating coded mask patterns based on Gaussian random fields (GRF) 
is proposed. In contrast to traditional algorithms based on cyclic difference sets, it is possible 
to construct mask patterns that encode a predefined point spread function (PSF). The viability 
of this approach and the reproducibility of the PSFs is examined, together with studies on the 
mean transparency, pixel-to-pixel variance and PSF deterioration due to partial shadowing. 
Sensitivity considerations suggest the construction of thresholded realisations of Gaussian 
random fields (TGRF) which were subjected to the same analyses. Special emphasis is given 
to ray-tracing simulations of the pattern's performance under finite photon statistics in the 
observation of point sources as well as of extended sources in comparison to random masks 
and the pattern employed in the wide field imager onboard BeppoSAX. A key result is that 
in contrast to traditional mask generation schemes, coded masks based on GRFs are able 
to identify extended sources at accessible photon statistics. Apart from simulating on-axis 
observations with varying levels of signal and background photon counts, partial shadowing 
of the mask pattern in the case of off-axis observations and the corresponding field-of-view is 
assessed. 

Key words: instrumentation: miscellaneous, methods: numerical, techniques: image process- 
ing 



1 INTRODUCTION 

In X-ray astronomy, focusing of radiation is so far feasible only 
for photon energies up to about 10 keV through grazing incidence 
reflection. Applied in Wolter-type mirrors, this method can pro- 
vide a very good angular resolution, i.e. down to C.'S in the case 
of Chandra^ and 4" - 12" for XMM-Newtorp- . The collecting area 
is maximised through the use of nested mirrors. The field-of-view 
(FOV) is limited by the grazing incidence condition set by the 
diffractive index of the mirror material to i 1°. At energies higher 
than 10 keV, focusing is technologically very hard to archieve. A 
workaround are coded mask imagers, where a position sensitive de- 
tector records the shadow of a mask pattern cast by the sources un- 
der investigation. The arrangement of sources can be reconstructed 
by cross-correlating the recorded shadowgram with the mask pat- 
tern. 

Coded masks have by now found a widespread use in high en- 
ergy astrophysics and there is a large number of successful missions 



such as BeppoSAX^ , currently flying intruments like INTEGRAL!^ 
and HETE-2^, and ambitious future projects, for instance SWIFJ*. 

In this paper, I propose coded mask patterns based on Gaus- 
sian random fields, because they enable the construction of a coded 
mask device for predefined imaging characteristics, i.e. for a given 
PSF. The shape of the PSF can be tuned to match the anticipated 
source profile. A beautiful example of a naturally occurring Gaus- 
sian random field is the pattern of fluctuations in the cosmic mi- 
crowave background (CMB). Analyses of WMAP data carried out 
among others by Cavon et al. ( 2001) and Komatsu et al. ( 2003) find 
the CMB consistent with Gaussian primordial fluctuations and have 
set upper limits on non-Gaussianity. 

After a recapitulation of coded mask imaging and existing 
mask pattern generation schemes in Sect.|5| GRFs are introduced 
in Sect.|3| The feasibility of GRFs in coded mask imagers is exam- 
ined in Sect.|4|with special emphasis on the performance of GRFs 
in realistic scenarios, i.e. under finite photon statistics and in the 
observation of extended sources (Sect. |5j. A summary of the key 
results in Sect.|5|concludes the paper. 
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2 CODED MASK IMAGING 

Coded mask cameras observe a source by recording the shadow 
cast by the mask onto the detector. The mask pattern is described by 
the position dependent transparency ifi(x). A shifted shadowgram 
- d tan(6)) is observed if the radiation incides under an angle 
6 with respect to the optical axis. The distance between the coded 
mask and the detector is denoted by d. The correlation function 
a(x), defined as 



a(x) = i//(x - d tan 6) i/)(x) = {if/(x - d tan + A)if/(A)},i, 



(1) 



peaks at xo = dland, from which the angle of incidence 9 = 
arctan(xo/(i) can be inferred. The PSF c(x), defined as the correla- 
tion function at normal incidence (6 = 0°), i.e. the auto-correlation 
function, reads 



c(x) = ij/(x) ij/(x) = {i//(x + A)i/^(A))a. 



(2) 



The influence of imperfections of the detector can be modelled by 
convolution of tf/(x) with suitable kernels describing the positional 
detector response (see, e.g. Schafer & Kawai 20^. Techniques for 
analys ing coded mask data have been summarised bv lSkinner et alJ 
h987h and Caroli et al. (1987) 

Random mask patterns as used in the HETE-2 satel- 
lite lin 't Zand et alj .1994.) consist of white noise. They are 
not ideal imagers, because their auto-correlation possess side- 
lobes and are not perfectly flat. Aiming at <5-like PSFs, 
mask patterns based on cyclic difference set s have been in- 
troduced by Gunson & Polvc hronopulo j Jl976 ^. As pointed out 
bv Fenimore & Cannon 1 1978) , these uniformly redundant arrays 
(URA) provide even sampling at all spatial scales. URA patterns 
are less susceptible to noise compared to truly random arrays and 
their auto-correlation function is a (5-spike with perfectly flat side- 
lobes in case of complete imaging. In this paper, I propose a 
method for constructing coded mask pattern encoding arbitrary 
PSFs. While the traditional masks are optimised for the observation 
of point sources, the PSFs of masks based on GRFs can be adjusted 
to the source profile of extended sources and make the observation 
of extended sources such as extended structures in the Milky Way 
possible. 



3 GAUSSIAN RANDOM FIELDS 
3.1 Definitions 

The statistical properties of a GRF are homogeneous and isotropic 
and the phases of diff^erent Fourier modes are mutually uncorre- 
lated and random. A consequence of the central limit theorem is 
then that the amplitudes follow a Gaussian distribution. Due to all 
correlations above the two-point level being either vanishing in the 
case of odd moments or being expressible in terms of two-point 
functions for even moments, the statistics of amplitude fluctuations 
in a GRF is completely described by its power spectrum P{k) (see 
eqn. 

Because the imaging characteristics of coded mask imagers 
are described by the PSF, which is defined to be the auto-correlation 
function of their mask pattern, i.e. by their power spectrum in case 
of isotropic PSFs, GRFs provide a tool for generating mask pat- 
terns with predefined imaging characteristics. The theory of struc- 
ture formation in cosmology and the description of the co smic mi- 
crowave b ackground makes extensive use of GRFs (c.f. IPeacockl 
Il999t lEonaair 1998). Their application is commonplace in gen- 
erating initial conditions for simulations of cosmic structure for- 



mation and in constructing mock CMB fields for simulating sub- 
millimetric observations. 



3.2 Algorithm 

Starting from the PSF c(jc), the Fourier transform C(k) is derived: 

C(k) = T[cix)]= f ^c(x)exp{-ikx). (3) 

The power spectrum P(k) is defined as the Fourier-transform C(k) 
of the auto-correlation function c(x). In more than one dimen- 
sion, an average of the Fourier transform C(k) of the statistically 
isotropic random field c(x) over all directions of the wave vector k 
at fixed length k = \k\ needs to be performed: 



P(k) = <|C(A:)|)|,|=,. 



(4) 



All elementary waves exp(ikx) with wave vectors in the k- 
space shell [\k\ , \k + Ak\] contribute to the variance cr^ = P(k) 
required by the power spectrum on scale k = |^|. In discretis- 
ing, the amplitudes ^{k) are set such that their quadratic sum 
2<re[|*i i;t+A*|] ll'C^)!" matches cr^ with the only exception 'Vik = 0), 
which is set to zero in order to ensure a vanishing expectation value 
of the realisation i^{x). The normal modes exp(iA:jt:) are modified by 
a phase factor e.xp{2niq), where q e [0, 1) is a uniformly distributed 
random number. By inverse Fourier transform, the normal modes 
l'(A) are brought to interference which finally results in the reali- 
sation, the real part of which is denoted by if/(x): 



(5) 



i//(x) = %{T'^ VVik)]) = ^ d^'^ ^(A:) exp(ikx + 2mq)j 

Alternatively, one may require the additional symmetry 
l'(-A) = ^"(A) in Fourier space (the complex conjugation is de- 
noted by the asterisk), which forces the realisation to be purely real. 
The flow chart eqn. jSJ summarises all steps: 



C(k) 



(■■■)\k\=k 



P(k) 



• exp(27ri(j) 



4'(A:) 



(6) 



c{x) \p{x). 

Due to the periodic boundary conditions imposed by the 
Fourier transform, the resulting realisations of the Gaussian ran- 
dom field have cyclic boundaries, which is a desirable feature 
for coded mask patterns. For reasons of numerical accuracy, it is 
strongly recommended to use shells in A-space with varying thick- 
ness Ak oc \k\^\ such that approximately the same number of dis- 
cretised modes contributes to the variance required by the power 
spectrum P(k). 

3.3 Choice of the PSF 

Although the algorithm outlined in Sect. I3.2l is capable of gener- 
ating random fields t//(x) encoding any isotropic PSF c(x), PSFs 
should be shaped like Lorenzian functions ci(x) or Gaussian func- 
tions Cg(x). The parameter <Tj, describes the spatial extent: 



Cl(x) 



ca(.x) 



exp - 



2<rl 



(7) 



(8) 



The normalisation has been chosen such that the maximum 
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correlation strength at x = is set to one. In the realisation 
the variable cr^^, th at parameterises the PSF ca n be interpreted as a 
correlation length. Skinner & Grindlav (1993') have pursued a re- 
lated idea and have suggested coded masks with two spatial scales. 
In contrast, the realisations considered here have an entire spectrum 
of length scales. 



3.4 Scaling applied to the Gaussian random fields 

If one aims at employing GRFs in coded mask imagers, the field 
has to be scaled such that it assumes values ranging from 4i{x) = 
(opaqueness) to full transparency (i^(x) = 1). This scaling ensures 
that the full dynamical range between is used and the modulation 
of the shadowgram as strong as possible. Hence, the sensitivity is 
maximised. One could think of two different linear transformations, 
the most intuitive being: 



if/(x) - min {tf>(x)] 



(9) 



max {ij/(x)] - min {tl/{x)] 

With the symmetry condition max {if/(x)] = -min {tff(x)] being 
fulfilled, the mean transparency (iff* (x)) is equal to 1/2: The mean 
(if/ (x)) = vanishes by construction, because each normal mode 
cos(kx) has a vanishing expectation value. In general, the realisa- 
tion tf/(x) will not fulfill the above mentioned symmetry condition. 

Instead, the scaling 



if>(x) if>'(x) : 



max(|i/r(x)|) 



-I- 1 



(10) 



ensures {if/') = 1/2 and will be used in the remainder of the pa- 
per. It should be noted that none of the above scalings strictly con- 
serves Gaussianity, because each particular realisation is scaled by 
its maximal amplitude and consequently, high amplitudes do not 
appear any more in an ensemble of realisations. 

Now that the mean transparency (i//') is fixed, the absolute flux 
from a source can be inferred from the number of measured pho- 
tons. The scaling eqn. may be taken advantage of in designing 
a mask that blocks a larger or smaller fraction of photons than the 
generic fraction of 1/2: In anticipation of Sect. l4.3l in the case of a 
realisation of a GRF encoding a Gaussian PSF ca(x) with tr,. = 8 
pixels, the probability density p{{tf/*)) d{i//*) of the mean trans- 
parency t = (if/*) is described by a Gaussian distribution with mean 
^i, = 0.504±0.082 and standard deviation cr, = 0.028±0.006 at 95% 
confidence. When constructing realisations of Gaussian fields for 
coded mask instruments, one obtains patterns with transparencies 
(i//*) e [ju, - cr,,//, + cr,] with a probability of erf(I/ V2) ^ 0.6827. 



3.5 Gaussian random fields for circular apertures 

For coded-mask experiments with a circular aperture it is possible 
to construct GRFs with azimuthal symmetry, in the same way as 
hexagonal uniformaly redundant arrays (H URA) are an adaptatio n 
of the URA patterns to circular apertures jFinser & Princ3ll985l) . 
Instead of constructing a GRF with plane waves as the solutions 
of Laplace's equation Aif/(x) = (fl^ -I- d^^ tf/(x) = in Cartesian 
coordinates (x = (x, y)) with boudnary conditions if/(x = -L) = 
= il/{x = L) (2L denotes the pattern's side length) one would 
resort to solving Aif/(r) = (fl; -I- l/rd,- + llr^d^^ipir) = in polar 
coordinates (r = (r, (j))) with the boundary condition i^(r = i?) = 
V0, where the radius of the aperture is denoted as R. i// is easily 
found as the solution to Bessel's differential equation and reads as: 





^Cm(r,(p) = J,„ [r ■ Z„,(e)] ■ exp(im(l>). 



(11) 



Figure 1. Normal modes il/f,„(r, </>) used for constructing GRFs for circular 
apertures, for (f, m) = (2, 3) (left panel) and (C, m) = (3, 4) (right panel). 



where the numbers { and m are only allowed to assume integer 
values. Z„,(f) is the f'" zero of the Bessel function J„. In Fig.Q 
two solutions are depicted for (I, m) = (2, 3) and ({, m) = (3, 4). 
In reality, it might be cumbersome to construct a GRF on the basis 
of the normal modes given by eqn. <11> due to Bessel function's 
complicated orthonormality relations. 



4 RESULTS 

In order to provide a visual impression, two GRFs encoding the 
above stated PSF with their auto-correlation functions are presented 
(Sect. 14. 1> . Subsequently, the reproducibility of the chosen PSF 
(Sect. 1^5), the pixel-to-pixel variance (Sect. R31 . the Gaussianity 
of the distribution of pixel amplitudes (Sect. l4!4l and the shape of 
the PSF under partial shadowing (Sect. l431 are examined. Finally, 
thresholded GRFs are introduced and the deterioration of the PSF 
of such thresholded realisations (Sect. l4.6l is addressed. 



4.1 Visual impression 

Following the above prescription, 100 realisations of GRFs encod- 
ing Gaussian and Lorenzian PSFs of different widths cr^ were gen- 
erated on a 2-dimensional square grid with 256^ mesh cells. Figs.|5| 
and |3| show a realisation of the GRF and its auto-correlation func- 
tion for a Gaussian and a Lorenzian PSF, respectively. In order to 
facilitate comparison, the widths of the PSFs have been chosen to 
be the same: cr^ = S pixels. The random fields are scaled to mean 
values of 1/2 (by means of eqn. <10» and the central correlation 
strength in the auto-correlation functions is equal to 1. The con- 
tours have a linear spacing of 0.1. The auto-correlation functions 
have the symmetry property that tl/(x) ifi(x) = ipi-x) ® ifi(-x). In 
the derivation of auto-correlation and cross-correlation functions, 
the balanced correlation scheme was used. The correlation func- 
tions were derived for ideal detectors, i.e. finite position resolution 
or similar imperfections were neglected. 

In comparing the realisations in Figs.|5|and|^one notices the 
larger abundance of small scale structures in the realisation encod- 
ing the Lorenzian PSF Cl(x) in comparision to the realisation de- 
rived for the Gaussain PSF cq(x)- This can be explained by the fact 
that the power spectrum Piik) declines oc exp(-k) and thus much 
slower than the power spectrum Paik) oc exp(-A;'). Both realisa- 
tions have been derived with the same random seed, i.e. the relative 
phases are identical and one immediately recognises similar struc- 
tures in t//a(x) and ifiL(x). 
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Figure 2. A realisation of a GRF i/zai^) (upper panel) for the Gaussian PSF 
cq(.x) and the auto-correlation function ipaW ® 'I'Gi^) (lower panel). 



4.2 Reproducibility of the PSF 

An important issue is the reproducibility of a chosen PSF c(x) 
in realisations generated with differing random seeds. This can 
be assessed by determining the auto-correlations of the scaled 
GRFs if/'(x) for all realisations within the data sample. In Fig.|4| 
the Gaussian target PSF co(x) and the auto-correlation functions 
ili'f,{x)®tlf'g{x) following from two realisations 4i{x) are shown. The 
error bars denote the sample variance derived from 100 realisa- 
tions of if/aix) following from different random seeds. The width 
of the PSF was chosen as cr^ = 8 V2 pixels for better visibility. 
Fig.|5]shows the analogous for the Lorenzian target PSF Cl(x) with 
CT;, = 8 V2 pixels. 

As Figs.|31and|5|illustrate, the functional shape of the target 
PSF c(x) can be reproduced with high reliability and the ratio of 
the peak-height to the correlation noise is ^ 40. However, there 
are minor imaging artefacts, namely very weak sidelobes: This is 
readily explained by the fact that the Fourier transform of a well 
localised PSF in real space is extended and affected by the cutoff at 
the Nyquist frequency ^Nyqiiist. which induces a sin(fcNy^„isf'^)/-«^-like 
modulation. Consequently, the sidelobes are suppressed in PSFs 
with large 0"^. The Lorenzian PSF is a bad choice in comparison 
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Figure 3. A realisation of a GRF iptM (upper panel) for the Lorenzian PSF 
Ci(,r) and the auto-coiTelation function iitix) (lower panel). 
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Figure 4. Cross section along the .t-axis through the central part of the 
auto-coiTelation function i/<q(x)8ii/<q(j:) for two different realisations 4i'q(x) 
(dashed) and the Gaussian target curve cq (x) (solid). 
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Figure 5. Cross section along the A:-axis through the central part of the 
auto-correlation function i//'^(x) if/'^(x) for two different realisations if/'^(x) 
(dashed) and the Lorenzian target curve ct(x) (solid). 



to the Gaussian PSF, because its Fourier transform C(k) oc exp(-A') 
decays slower and is consequently more affected by the cutoff at 
^Nyquist- Interpreting cr,- as the correlation length of the GRF, it is 
clear that in the limit of very narrow PSFs cTv assumes very small 
values, i.e. the amplitudes tf/{x) for neighbouring pixels start loos- 
ing their correlation. This, however, does not correspond to white 
noise masks because the amplitude distribution is still Gaussian 
(c.f. Sect. l4.4> and not bimodal, as in the case of white noise masks. 

Due to the high confidence with which a chosen PSF is repro- 
duced, the number of realisations to be examined is very small. On 
the contrary, relying on truly random patterns, the number of neces- 
sary realisations with the accompanying tests may be very high: For 
HETE-2, where such a random pattern is used, 10' realisations had 
to b e generated that were subjected to certain boundary conditions 
(see lin 't Zand et alll994l) . 



4.3 Pixel-to-pixel variance 

In sensitivity considerations carried out bv lin 't Zand et al] il994) 
for purely random masks, i.e. masks consisting of either transpar- 
ent (tff'(x) = 1) or opaque (if''(x) = 0) pixels, optimised mean trans- 
parency (i^'} and standard deviation - {ij/'Y are derived to 
be equal to 1/2. In that way, the variance and therefore the modula- 
tion of the signal is maximised. For the GRFs considered here, the 
variance and hence the modulation of the shadowgram is notica- 
bly smaller. In TableQ the mean transparencies the variance 
{ijj'^) - (lA')^ ™d the standard deviation — (i//')- together 
with their respective uncertainties for a set of GRFs encoding Gaus- 
sian PSFs with differing width cr,- are summarised. 

One would expect that with increasing PSF width cr^ the vari- 
ance decreases, which would be explained by the fact that the vari- 
ance is given by a weighted integration over the power spectrum 
P(k). For increased position resolution, i.e. a narrow PSF c(x), a 
wide power spectrum P(k) is needed, which in turn would lead to a 
high variance. 

This simple argument however, does not straightforwardly ap- 
ply to the scaled realisations at hand: As laid down in eqn. <10> . 
the field if/{x) is modified by a factor depending on the maxi- 
mal value \tf/(x)\ of the particular realisation. The occurence of a 



PSF width 


mean transparency 


variance 


standard deviation 






<^'2) - <^')2 


V<iA'2) - <^')2 




1/2 


0.013 ± 0.002 


0.116 ±0.009 


4 


1/2 


0.014 ± 0.002 


0.121 ±0.010 


4V2 


1/2 


0.016 ±0.003 


0.128 ±0.012 


8 


1/2 


0.018 ± 0.003 


0.135 ±0.013 



Table 1. The mean transparencies (lA'), the variance {1//'^) - (1//')^ and the 
standard deviation ^j(lf^'-) - (if/')^ together with their respective uncertain- 
ties (Icr) for a set of GRFs encoding Gaussian PSFs with differing width 
o-v. 



high amplitude is following a Gaussian distribution with variance 
oc f d-k P(k). This means, that in the case of narrow PSFs c(x), i.e. 
for wide power spectra P(k) , the field il/(x) is more likely to assum e 
large amplitudes (compare Cartwright & Longuet-Hi esins 1956). 
The latter effect is of great importance and causes the surprising 
result that the measured variances in if/'(x) are larger for extended 
PSFs. 

Comparing coded masks based on GRFs with purely random 
fields, the modulation of the shadowgram decreases by a factor 
~ 3 ... 4. Therefore, the sensitivity is expected to be weaker. While 
the above consideration is only valid for the observation of point 
sources, sensitivity is most likely to be gained in the observation 
of extended sources. For those sources, it is possible to adjust the 
PSF to the expected source intensity profile. In this case, modula- 
tions below the scale of the object to be observed are discarded - 
this corresponds to applying Wiener filtering to the recorded shad- 
owgram prior to source reconstruction. 

4.4 Distribution of the pixel amplitudes 

As Fig.|6|illustrates, the pixel amplitudes i/'(x) follow a Gaussian 
distribution, irrespective of the encoded PSF, 



p(iA')diA' 



^J2no 



■ exp 



(12) 



as a co nsequence of the central limit theorem (see lKendall & Stuart! 
il958l» . The mean and variance of that particular realisation have 
been determined to be ji^ = 0.5000 ± 0.001 and cr^ = 0.1277 ± 
0.0007 at 95% confidence. For illustrative purposes, a Gaussian 
PSF with cr, = 2 V2 pixels has been chosen. 

Again, it should be emphasised that the scaling eqn. )10t . 
while being reasonable from the physical point of view, is not con- 
serving Gaussianity. This is for the application at hand not a seri- 
ous limitation, because the variance of the distribution p (tf/') dtf/' is 
small compared to 1 . 



4.5 Partial shadowing 

It is interesting to see how partial shadowing affects shape and am- 
plitude of the auto-correlation function. If a source is observed at 
large off-axis angles, the shadowgram cast by the coded mask onto 
the detector is incomplete and reconstruction artefacts emerge in 
the correlation function. In order to examine the extent to which 
the PSF suffers from partial shadowing, the amplitudes if/cix) in a 
margin amounting to a fraction of 25%, 50% and 75% of the to- 
tal area have been set to zero and the cross-correlation function 
t//o(x) if/-^'"'''"' (x) has been determined with the full coded mask. 
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Figure 6. Probability density p (if/') dif/' of the pixel amplitudes i//' (x) (cir- 
cles) and the best-fitting Gaussian for a particular realisation of a GRF. The 
error bars are Poissonian errors. 
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Figure 8. Reconstructed auto-correlation function ifi'g(x)^i//'^'"^"^(x) for a 
shadowgram of which only ^3% have been imaged onto the detector. 



1.4 



1.2 - 




_Q_2l 1 1 1 1 1 1 

-60 -40 -20 20 40 60 

pixel coordinate x 

Figure 7. Cross-correlation function i//'^(x)^i//-^'"'°''''' (x) and the respectively 
expected PSF c(x) with a shadowed margin corresponding to 25% (dashed, 
squares), 50% (dotted, diamonds) and 75% (dash-dotted, crosses) of the 
total area and, for comparison, the PSF for full imaging (solid, circles). 



As Fig. shows for a Gaussian PSF with cr,. = 8 pixels, the 
PSF drops in central amplitude according the unshadowed area, but 
otherwise its shape remains unaltered. A second observation is that 
the amplitude of the sidelobes is unaffected by the partial shadow- 
ing. 

The reconstructed PSF 4i'q{x) ® tf/*'"'°"Xx) for the case of ra- 
diation from a source situated at large angles away from the opti- 
cal axis, where only 1/32 of the mask has been imaged onto the 
detector is depicted in Fig.[8| Even though a tiny part of the masls 
amounting to ^ 3% has been imaged, the correlation peak is clearly 
recognisable and its peak value is a factor =^ 4 above the correlation 
noise. 



4.6 Thresholded realisations 

Due to possible technical complications in attempting to build a 
coded mask pattern based on a GRF with quasi-contiuous opaque- 
ness, thresholded realisations are considered. A second argument in 
favour of thresholded realisations would be their achromatic prop- 
erties, because the mask has to be constructed from the field if/'(x) 
for a specific photon distribution in order to assure the maximal 
modulation of the shadowgram cast onto the detector. Yet another 
argument in favour of thresholded realisations of GRFs is their bet- 
ter sensitivity, because they imprint a stronger modulation of the 
shadowgram compared to smoothly varying GRFs. 

In thresholded realisations, mask elements are taken to be 
transparent, if the value t//(x) of the realisation is greater than zero, 
conversely, for values if/(x) < the mask element is set to be 
opaque. An example for a thresholded realisation of a GRF and 
its PSF is given in Fig.|9| 

An important issue is the degradation of the PSF i^g (x) ® 
(jc) imposed by the thresholding. As Fig. llOl illustrates. the re- 
sulting auto-correlation function is pointy and its kurtosis is larger 
than zero (leptokurtic). This results from the fact that small scale 
power is added by the thresholding: In order to construct a step 
transition, more small-scale Fourier modes are needed, which leads 
to an additive power law contribution oc k^^ in the power spectrum 
P(k), such that the power spectrum acquires Lorenzian wings. The 
point spread function i//q(x) ® if/'^ix), being the inverse Fourier 
transform of P(k), can then be approximated by two decaying 
branches of an exponential, which readily explains the pointiness. 
The target Gaussian PSF Cg(x) with cr^ = 8 V2 pixels is shown 
for comparison. Again, the error of the auto-correlation function is 
estimated by determining the sample variance in 100 realisations. 

The size distribution of the patches as a function of threshold 
value can be described by means of the P ress-Schechter theory wel l 
known in cosmology, ^ress & Schechteil iT97l . lBond et ai]jl99ll) 
and lMo & White! i2002l) provide the mathematical foundation. 
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Figure 9. A thresholded realisation ^^q(x) of a GRF ^a{x) (upper panel) 
for the Gaussian PSF cg(jc) and the auto-correlation function i/r^ (a:)®!/'^ (jr) 
(lower panel). 



5 RAY-TRACING SIMULATIONS INCLUDING FINITE 
PHOTON STATISTICS 

Extensive ray-tracing simulations were performed describing the 
imaging of point sources with a finite number of photons 
(Sect. 15. 1> . and the attainable sensitivity in such an observation 
was assessed (Sect. 15. 2t . The analogous was carried out for the 
observation of extended sources (Sect.|53}. Finally, the size of the 
field-of-view in the case of GRFs compared to traditional masks is 
examined (Sect. l7!4l . In the following, coded mask patterns based 
on Gaussian random fields are compared to purely random mask 
patterns and the mask pattern used in the WFI-instrument onboard 
BeppoSAX. 

5.1 Simulation setup 

In the following, the performance of the coded mask is examined 
as a function of photon statistics. The statistical significance cr of a 
simulated observation is defined to be 



o 

3 0.4 




(13) 



Vbg 



pixel coordinate x 

Figure 10. Deterioration of the auto-correlation function iP'q(x) 8 if/'gix) 
from two thresholded realisations iI/q(x) (dashed) in comparison to the ini- 
tial Gaussian PSF ca(.x) (solid). 



where A'source and A'bg denote the source and background count 
rates, respectively. Here it should be emphasised, that cr, A^somcc and 
A'bg always refer to the number of actually detected photons which 
makes a difference when considering the coded mask employed in 
BeppoSAX's WFI instrument, in which the average transparency is 
not equal to 1/2. 

Observations were simulated by randomly choosing 2A'soiircc 
homogeneously distributed photon impact positions x across the 
mask face. In order to emulate the random process of photons pen- 
etrating the mask, a homogeneously distributed random number r 
from the interval re [0 ... 1] was drawn for each photon, and com- 
pared to the value ifi(x) of the GRF at the same position x. In the 
case r > tf/(x) the photon was assumed to be able to penetrate the 
mask, whereas in the case r < ifi(x) the photon was taken to be ab- 
sorbed by the mask. For BeppoSAX's mask pattern, which has an 
average transparency of 1/3, a total number of SNsoma: photons was 
simulated. 

For the background, which was assumed to be homogeneous, 
A'bg photon impact positions were determined and the count rates in 
the corresponding pixels were increased accordingly. Background 
count rates were fixed to a value of A'bg = lO'* photons, which are 
typical for an instrument like WFI in a 100 second exposure. 

The resulting field ifi'-"'"^ containing the number of photons 
that struck a certain pixel was then correlated with the origi- 
nal mask pattern if/, again using balanced correlation. In the next 
step, the highest peak was localised in the correlated data field 
^ significance E was determined by comparing the 

peak height a„,ax = maxji^*^""' ® if/] to the level of fluctuations 
0-2 = {{if/('""> (g) ) in the field. If the peak had a significance 
^ = o^a^/cr^ exceeding 3 and was located at a position which de- 
viated less than half a PSF width from the nominal position, the 
simulated detection was taken to be successful. A particular real- 
isation of a Gaussian random field was exposed to 100 simulated 
photon distributions from which the detection probability p (i.e. the 
occurence of a > Scr^-peak located at the correct position) and the 
false detection probability q (i.e. the occurence of a > 3cr,;,-peak at 
a wrong position) was derived. The sample variance in comparing 
100 realisations of Gaussian random fields was used to derive er- 
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Figure 11. Point source sensitivity in on-axis observations of a set of 
GRFs: Tlie detection probability pia) (solid lines, closed symbols) and the 
false detection probability q(<T) (dashed lines, open symbols) are plotted 
as functions of statistical significance o" for PSF widths o"j = 2 (squares), 
iTv = 2 V2 (circles), cr, = 4 (triangles) and cr,- = 4 V2 (diamonds), in com- 
parison to purely random masks (dotted line, stars) and BeppoSAX-WVl 
pattern (dash-dotted line, crosses). In contrast to the ensemble of GRFs it is 
not possible to state an ensemble variance for p(o") and q{tT) in the case of 
BeppoSAX's pattern. The data points have been slightly displaced for better 
visibility. 

rors on p and q. For the purpose of this work, the detector efficiency 
and position response were assumed to be ideal. 

5.2 Point source sensitivity of a set of Gaussian random fields 

Fig. II ll shows the detection probabifity p and the false detection 
probability ^ as a function of photon statistics, expressed in terms 
of statistical significance cr for GRFs, a purely random mask and 
BeppoSAX's URA pattern. The source was assumed to lie on the 
optical axis, i.e. the mask pattern is imaged completely onto the 
detector. Common to all mask patterns is the fact that p rises with 
statistical significance, and that q drops accordingly. But while reli- 
able observations can be done using the BeppoSAX-pattem or ran- 
dom patterns even at low photon statistics of 2 - 3<t, the patterns 
based on GRFs require high photon fluxes. For them, observations 
are feasible starting from ^ 9tr. The reason why GRFs are less sen- 
sitive to the traditional mask pattern is the fact that they imprint a 
weaker modulation of the shadowgram. Furthermore, one immedi- 
ately notices the trend that the patterns are more sensitive for wider 
PSF widths due to the increase in variance of the mask pattern with 
increasing PSF width. Thus, position resolution is traded for sensi- 
tivity. 

Fig. 1121 shows the analogous results for an off-axis observa- 
tion in which only half of the mask pattern has been imaged onto 
the detector. The result corresponds to the findings for the case of 
normal incidence, but p(cr) and q(cr) are shifted to higher values of 
cr, which is due to the fact, that only half of the photons actually 
reach the detector and that the reconstruction has to cope with the 
decreased signal. Again, one attains higher sensitivities for wider 
PSFs in the case of patterns based on GRFs. 

Common to all figures is the fact, that the curves p(cr) and g(cr) 
are not adding up to one, which is caused by the combined criterion 
where apart from the correct peak position a minimal peak height 
above the correlation background is required, which is often not 
fulfilled in the cases of low photon statistics. 
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Figure 12. Point source sensitivity in off-axis observations (50% obscura- 
tion) of a set of GRFs: The detection probability pier) (solid lines, closed 
symbols) and the false detection probability q((T) (dashed lines, open sym- 
bols) are given for PSF widths o", = 2 (squares), cr^ = 2 V2 (circles), 
(Tj- = 4 (triangles) and ctj = 4 V2 (diamonds), in comparison to purely ran- 
dom masks (dotted line, stars) and BeppoSAX's WFI pattern (dash-dotted 
line, crosses). 

5.3 Sensitivity in observations of extended sources 

In addition, suitable simulations were carried out in order to assess 
the performance of GRFs in the observation of extended sources, 
such as supernova renmants, structures in the Milky Way and clus- 
ters of galaxies. Typical sizes of those sources range between ar- 
cminutes and a degree. For simplicity, the source was assumed to 
be described by a Gaussian profile with extension o-pyame = 2 pix- 
els. The shadowgram recorded in observations of extended sources 
are superpositions of slightly displaced point source shadowgrams, 
where the relative intensities follow from the source profile. Conse- 
quently, the imaging of extended sources is simulated by convolv- 
ing the mask pattern with the source profile prior to the ray-tracing. 
Despite that, the image reconstruction has been carried out with the 
unconvolved mask pattern. 

Fig. 1131 gives the dependence of the detection probability p 
and the corresponding q on the photon counting statistic cr. In the 
observation of extended sources, the patterns based on GRFs are 
superior to the traditional approaches: While reliable detections can 
be achieved starting from cr 10 (for cr^ = 4 V2) up to cr * 20 
(for cr,. = 2), the performance of the traditional masks is notably 
worse. At the examined levels of photons statistics, the detection 
probability p stays close to zero and shows but a shallow increase 
with cr in the case of BeppoSAX's URA pattern. 

The good performance of the GRFs, and their decreasing per- 
formance with correlation length, i.e. PSF width cr, is of course 
to be traced back to the fact, that mask patterns with large struc- 
tures are less affected by the convolution with the source profile 
than mask patterns exhibiting small structures; in the extreme case 
of random masks or BeppoSAX's pattern, the structures are washed 
out and consequently, the modulation of the shadowgram is very 
weak. This can be circumvent, however, by tuning the angular size 
of a mask pixel to match the angular size of the source to be ob- 
served. 

5.4 Field-of-view in tlie observation of point sources 

Now, the size of the field-of-view, i.e. the minimal fraction of the 
mask pattern required to be imaged onto the detector in order to 
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Figure 13. Sensitivity in on-axis obsei^ations of extended sources of a set 
of GRFs: The detection probability p(<T) (solid lines, closed symbols) are 
given along the false detection probability q{(T) (dashed lines, open sym- 
bols) for PSF widths of cr, = 2 (squares), ctj- = 2 V2 (circles), tr^ = 4 
(triangles) and cr, = 4V2 (diamonds), in comparison to purely random 
masks (dotted line, stars) and BeppoSAX's WFI pattern (dash-dotted line, 
crosses). 
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Figure 14. Decrease of the field-of-view: The source detection probability 
p(Q) (solid lines, closed symbols) and the false detection probability q(Q) 
(dashed lines, open symbols) as functions of the obscuration Q are shown 
for a set of Gaussian random fields for varying PSF width: <t^ = 2 (squares), 
o",. = 2 V2 (circles), 0", = 4 (triangles) and CTy^ = A V2 (diamonds). In 
comparison, a purely random mask (dotted line, stars) and BeppoSAX's WFI 
pattern (dash-dotted line, crosses) are considered. 



yield a significant detection peak is investigated. For that purpose, 
the point source detection probability p and the false detection 
probability q are considered to be functions of the obscuration Q, 
which is defined as the fraction of the mask area imaged onto the 
detector. The number of background photons was kept fixed to be 
A^bg = 10"*, while the number of source photons A'somce was dimin- 
ished by this factor of Q prior to the ray-tracing. Their number was 
fixed to yield a significance of cr = 20 for 2=1, i.e. for the case 
of complete imaging. The background photons were assumed to be 
homogeneously distributed. The simulation and the derivation of 
the values for p{Q) and q(Q) were carried out in complete analogy 
to Sect.lO 

The results are depicted in Fig. 1141 While the traditional pat- 
terns show a good performance and have a high detection probabil- 
ity piQ) for values of g i 0.1 (BeppoSAX's pattern) and Q i 0.2 
(random mask), the GRFs fall behind significantly in performance. 
Imaging is only possible in the cases where a fraction of at least 
Q = 0.5 . . . 0.6 of the mask has been imaged onto the detector, re- 
sulting in a decrease of the field-of-view of about a factor of 3 ... 5, 
which renders the usage of GRFs very unlikely in survey missions. 
Again, the GRF patterns encoding wide PSFs are more sensitive 
and yield larger fields-of-view than GRFs with narrow PSFs. 



6 SUMMARY AND OUTLOOK 

In this article, a new algorithm for generating coded masks is pre- 
sented that allows the construction of a mask with defined imaging 
properties, i.e. point spread functions. 

• The viability of constructing a coded mask for a predefined 
PSF as a realisation of a GRF has been shown. For realisations 
generated with differing random seeds, the shape of the PSF is re- 
producible with high accuracy. Due to the reproducibility of the 
PSF, the parameter space is greatly reduced and the necessity of 
running extensive Monte-Carlo simulations is alleviated. 

• The generation of 2-dimensional URA patterns requires the 
number of pixels in each direction to be incommensurable, i.e. they 
are not allowed to have a common divisor While twin prime num- 



bers exist, mask patterns generated that way are almost, but not 
quite square l.Mivamoto..l977; .Proctor et al. 1979) . Coded masks 
based on GRFs may have any side length and any ratio of side 
lengths. Additionally, sizes chosen equal to 2", « 6 N enable the 
usage of fast Fourier transforms. Realisations of GRFs have cyclic 
boundary conditions which is a desirable feature for coded mask 
imagers. 

• The average transparency of coded mask patterns based on 
scaled GRFs is equal to 1/2, irrespective of the PSF they encode. 
The pixel amplitudes of a realisation are Gaussianly distributed as 
a consequence of the central limit theorem. The pixel-to-pixel vari- 
ance, however, is smaller in the case of GRFs compared to purely 
random fields, which results in a weaker modulation of the shad- 
owgram and hence the sensitivity is expected to be smaller. The 
variance shows the trend of decreasing with increasing PSF width, 
which is caused by the scaling with the maximal values of the real- 
isation. 

• Coded masks based on GRFs are chromatic in contrast to 
purely random fields: The mask pattern has to be designed for a spe- 
cific spectral distribution of photons due to semi-transparent mask 
elements. Any mismatch in the photon spectrum of a source un- 
der observation would result in a less pronounced modulation of 
the shadowgram, which in turn affects the sensitivity of the coded 
mask imager. A possible workaround is the usage of thresholded 
Gaussian random fields, that show pointy auto-correlation func- 
tions in contrast to smooth target PSFs. Another advantage is their 
enhanced sensitivity due to the stronger modulation of the shadow- 
gram. The properties of thresholded realisations, however, show a 
large sample variance which requires selections with suitable crite- 
ria after construction. 

• Ray-tracing simulations including finite photons statistics and 
background noise show, that the sensitivity of GRFs falls behind 
that of purely random masks and URA patterns like the one em- 
ployed in BeppoSAX by a factor of 2 ... 3 in the observation of 
point sources, depending on PSF width. For GRFs, the sensitiv- 
ity was found to depend exponentially on PSF width, one is trading 
sensitivity for position resolution. 
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• The sensitivity of patterns based on GRFs is significantly bet- 
ter in the observation of extended sources because their comparably 
large structures are less affected by the convolution with the source 
profile than traditional masks that possess pronounced structures on 
small scales. 

• Finally, the size of the field-of-view of GRFs in comparison to 
traditional masks is examined. It is found that reliable imaging can 
only performed with GRFs, if a large fraction of the mask is imaged 
onto the detector. In contrast, purely random masks and especially 
BeppoSAX's URA pattern enable imaging at large off-axis angles. 
Comparing the resulting fields-of-view for the preset number of 
photons shows, that the field-of-view of patterns based on GRFs 
are smaller by a factor of 3 ... 5 (depending on PSF width). 

Although the shortcomings of Gaussian random fields with 
respect to point source sensitivity, chromaticity and localisation ac- 
curacy make their usage in observing point sources doubtful, they 
may find application in observations of extended sources, while si- 
multaneously providing a moderate performance in the observation 
of point sources. Coded mask patterns on the basis of GRFs may be 
aesthetically pleasing because they utilise an abstract cosmological 
concept for a technological appUcation. 
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